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Abstract. We evaluate in closed form, for the first time, certain classes of double series, which are remindful 
of lattice sums. Elliptic functions, singular moduli, class invariants, and the Rogers-Ramanujan continued 
fraction play central roles in our evaluations. 



1. Introduction 

In this paper we establish elementary evaluations of certain 2-dimensional infinite series. 
For example, 



oo oo 



n=— oo m=— oo 



(5m)2 + (5n+ 1)2 5v^ V 

+ ^log(ll + 5v^' 



y y ^l'"^" ^_^iogfv^ + i-v^5 + 2v^ 



(1.1) 



which is a problem submitted to the American Mathematical Monthly [|T5l . The algebraic 
numbers on the right-hand side of (11.11) arise from special values of the Rogers-Ramanujan 
continued fraction. In general, elementary evaluations are quite rare for higher-dimensional 
lattice-type sums. For instance, the third author has examined both double and quadruple 
sums in connection with Mahler measures of elliptic curves; those sums typically reduce to 
values of hypergeometric functions ifTOl . [[TT|. [|T2|. The most famous higher-dimensional 
sum is the Madelung constant from crystallography BH, 0, jH, [[HI. It is highly unlikely 
that Madelung's constant possesses an evaluation in closed form. 

We produce many additional results along the lines of (|l.ll) . In fact, we show that it is 
possible to evaluate 

Fia,b){x):= J2 (xm)2 + [an + hf ^^'^^ 



n=— oo m=— oo 



for any positive rational value of x, and for many values of (a, 6) G N^. Since the se- 
ries is not absolutely convergent, we will calculate the n-index of summation using = 
lim^v^oo X]-Af<n<Af When a < 6, results from Ramanujan's notebooks apply, and the val- 
ues of the sums can be deduced from classical results in theta functions and g-series. When 
a > 6, the situation is slightly more complicated. In those cases we require the added hy- 
pothesis that a E 2Z, and then we use properties of Jacobian elliptic functions. 
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2. The main theorem 

Before proving our main theorem, we note that it is possible to calculate -F(a,b) (x) to high 
numerical precision with the formula 

^ ^ f-D^-csch 



7r(a?i+6) 

FiaM^) = - E — TT^' (2.1) 

X ^-^ an + 



n=—oo 



which can be established by substituting the partial fractions decomposition for csc\i{z) BH 
p. 28, Entry 1.217] in (11.21 ). Formula (12.11) provides an easy way to numerically verify results 
such as (11.11 ) and (13.201 ). For instance, we calculated F(5_i)(l) to more than 100 decimal 
places by summing over — 15<n<15. 

Theorem 2.1. Suppose that a and b are integers with a > 2, (a, b) = 1, and assume that 
Re(x) > 0. Then 

„ a-l / oo \ 

F^aM^) = -- 5^a;-(2^+i)Mog n (1 - ^''^y"'^') (1 - uJ-'^-V'''^') , (2.2) 

j=0 \m=0 / 

where u = e'^*/" and q = e~'^l^. 

Proof. Suppose that is a positive integer. Then using the transformation formula for the 
theta function ip{—q) = J2n=-ii~^)'^^"^' t^' P- Entry 27(ii)] and inverting the order of 
summation and integration twice by absolute convergence, we find that 

/ 1 \m+n roo I \ 

y 7 hr^, 7TT7 y (-1)" / e-"('^"+'')'" ,^yn^-.m^x^u 

^ (xmy + (an + b)^ ^ 7o lA^^ ^ i 

meZ ^ J ' \ 'J -N<n<N ^' VmeZ / 

-N<n<N 



=7r 



-N<n<N \^y'^m=0 



-iY. (-1)" 



e <<'^+^) «e — 



X — ' — ' /n \ U 

m=0-N<n<N 

The substitution of a standard K-Bessel function integral ^ p. 384, formula 3.471, no. 9] 
leads to 



Ul 



(-1)-+- _27ry. (_l)n^(2m+l)|an+6| 

^ (xmY + (an + bY x ^ Ian + 61 

meZ ^ I ' \ 'I m=0 -N<n<N ' ' 

-N<n<N 

Notice that the condition Re(s) > ensures convergence. Now recall that to = e'^*/". If 
a > 2, (a, b) = 1, and |r| < 1, then 

n=— oo ' ' i=0 
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which is easily verified by comparing Taylor series coefficients in r. Notice that (12.41 ) in- 
volves an infinite series in n, whereas (12.31 ) imposes the restriction that n E {—N, N). If 
we divide the left-hand side of (12.41) into three components, — oo<n<— A^, — A^<n< 
N, N < n < oo, and use a crude error estimate to bound the terms where n > N and 
n < —N, then we find that 

S l + = log ((1 - .''-'r) (1 - c.-<=i«)r)) 

-N<n<N ' ' j=0 (2 5) 

+ ' 



Substituting (12.51) into (12.31) leads to equation (12.21) plus an error term. The error term can 
easily be seen to approach zero as — )■ oo. □ 

3. SIMPLIFICATION FOR a G {3, 4, 5, 6} 

Although (12.21) was not difficult to prove, the deduction of results such as (|l.ll) from (12.21) 
is usually more difficult. In this section we examine the cases where a E {3, 4, 5, 6}. In these 
instances we can assume that 6=1 without loss of generality, because a standard symmetry 
(e.g. n — 7- —n) can be used to recover the other possible values of F(^a,b)ix). The same 
symmetry immediately implies that F(^2,i){x) = 0. 

Let us briefly recall the g-series notation 

oo 
j=0 

Following Ramanujan's notation for theta functions, define 

oo oo 
n=— oo n=0 

xiq) = ?^)oo, /(-g) = (g; g)oo. 

We need the famous Jacobi triple product identity [3, p. 35, Entry 19] for Ramanujan's 
general theta function /(a, 6), and they are given by 

oo 

f{a,b):= ^ a^^'^+^^/^^^^^-i^/^ _ (-a;a6)oo(-&;a&)oo(a&;a&)oo, |a6| < 1. (3.1) 

n=— oo 

We note the easily proved evaluation [3, p. 34, Entry 18(ii)] 

/(-l,g) = 0, |g|<l. (3.2) 
The Rogers-Ramanujan continued fraction 

^(^)- — + Y + T + T + ---' I^K^' 

also plays an important role in this paper. We always assume that g^/^ takes the principal 
value, so that R{—q) assumes a real value if g G (0, 1). 
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Theorem 3.1. Suppose that q = e""/^. Letu = —q^^^xi.Q.)/x^{<l^)' " = 1— v?^(— g^)/v?'^(g*), 
and = R{—q)R^{q^). Then 

(3.3) 

^(4,1) (^) = - ^log f^^l , (3-4) 



V2x ''Vl + v^ 
F ix)- , / 2-/i + 18/i^ + /x3 + 2/x^ + 5V5(/i + /i^) 

5v^x \2 - /i + 18/i2 + + 2/^4 - 5\/5(/i + /i^) ' 

TT , / 1 + yU - yU^ 

log 



5x \1 — 4/i — 

p ^1 ^ ^(-g4)_3 ( 36)^2v/3g/(V^) ^ 

F(6iUa;) = ^log ^ . (3.6) 

^ v^x l^y.(-g4) - 3^(-g36) _ 2v^g/(-g24) ^ 

Proof. We begin by proving (13.31) . If we set (a, h) = (3, 1), then (12.21) immediately reduces 
to 

^(3.(.) ^ -I E COS (^) log n (l - 2 COS (^) . 

To finish the calculation, let US briefly take a = 1— (p^{—q) /(p'^{q) and /S = l—(p^{—q^)/(p'^{q^). 
Then by the inversion 

formula for g-^/24^(^^^) ^ p ^24, Entry 12(v)], we have 



It is known that a and /3 admit birational parametrizations a = p{2 + pY/ (1 + 2p)^ and 
/3 = p^{2 + p) /(I + 2p) [3, p. 230, Entry 5(vi)]. Thus we obtain 

f{l-p){2 + p)V'^ q^'\{q) fp{l+p)^'/' 



X\q) V 2(l + 2p)2 y ' x^(g3) 
Finally, it is easy to see that if u := — then 

/ 1 + w3 N V3 



,3 



Substituting this last result into (13. 7|) completes the proof of (13. 3|) . 
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Next we prove (13.41) . Notice that if (a, b) = (4, 1), then (12.21 ) becomes 

^<M.(-) —It - (^^) >og n (i - 2 cos (^(BLii) ) + 

j=0 ^ ^ m=0 V V / / 

Letting to = e'^^^^ and using the Jacobi triple product identity (|3.1I) . we find that the denomi- 
nator on the far right side of (13.81) is equal to 

oo 

F{q) := n (l + V^?'™^' + ?'"^') (1 - 9'™""') 

m=0 

= (-wg; g^)oo(-wg; g^)oo(g^; g^)oo 

^ n (3-9) 

= 2^ u q 

m=—oo 
oo 

n=— oo 

Because the initial infinite product is real-valued, the imaginary terms above sum to 0. Al- 
ternatively, this fact also follows from (|3.2I) . Hence, from (|3.9I) . 

^ oo 

* n=— oo 

oo 

n=0 

= V'l-g'') + V2#(-g')- (3.10) 

A similar argument provides a similar representation for the numerator on the far right side 
of (13.81) . Hence, using (13.101) and its aforementioned analogue in (13.81) . we are led to the 
closed form 

-rf4 1-)UC = ^ log ^ . 

If we take a = 1 — '^'^{—q^)/'^'^{q^) and 2; = ip'^{q^), this expression reduces to (|3.4I) after 
applying [H p. 122, Entry lO(iii)] and [H p. 123, Entry ll(ii)]. 

Next we prove (13.51) . This case is substantially more difficult than the previous two. If 
(a, b) = (5, 1), then (12.21) becomes 

j=0 ^ ^ m=0 V V / / 

= _iLioe^-^loe TT 1 - + g'"^ .3 

m odd 
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where /3 = ^-^^ and a = ■ The second equality in (13.1 II) was obtained by collecting 
terms of the form log(X). Now we use several entries from Ramanujan's lost notebook. 
By [2, pp. 21-22, Entry 1.4.1, Eqs. (1.4.3), (1.4.4)], 



m odd 



By |l2l p. 33], we can parameterize R^(—q) and R^{q^) in terms of = R{—q)B?{q^) with 
the identities 

(3.13) 

Therefore (|3.12l) becomes 



+ ^ / 2 - /i + 18/x2 + + 2/^4 + 5x/5(/i + ^3) ^ 



m odd 



Finally, notice that g) = (— g; g)oo- After replacing q by — g in [2, p. 37, Entry 1.8.5] 

and simplifying, we have 

X(g^) _ 1 + ^ - ^3 
(g) 1 - 4/i - ■ 

Substituting (13.141) and (13.151) into (13.111) concludes the proof of (1331) . 

The proof of (13.61) is similar to the proof of (13.41) . and we leave this calculation as an 
exercise for the reader. Note that the operative result 



n 

n=0 \ 



1 + y3g2"+l + g4n+2\ -^(-g^) + ?,<^{-q^^) + 2y/?,qf{-q 



24^1 



- V3g2"+i + g4«+2 / -^(-g4) + 3^(-g36) - 2v^g/(-g24) 



follows easily from the Jacobi triple product identity (13.11) . □ 

Now we derive some explicit examples from Theorem 13. 1[ All of our identities fol- 
low from well-known g-series evaluations. We begin with F(3_i)(a;). Notice by f2}, p. 95, 
Eq. 3.3.6], that u = G{—q), where G{q) denotes Ramanujan's cubic continued fraction 
defined by 

^, , g^/3 g + g2 q"^ + q^ q^ + q^ , , 

Gid) ■■= — + ^ + + + . . . ' 1^1 < 1- 

It follows that F(3^i)(x) can be evaluated by using formulas for G{—q). When x = 1 we 
appeal to ^ p. 100, Eq. (3.4.1)] to find that u = G(-e-^) = which yields 

oo CO , ■-m+n n„ , . 



n=—oo m=—oo 
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When X = ^ we appeal to E p. 101, Theorem 3.4.2]. We have u = G(-e-''^) = 
(x/5-3)(v^-73) ^ therefore 

E E n^ = ^logf8(4-v^)). (3.17) 



n=— OO m=— OO 



It is possible to obtain many additional formulas for F(3_i)(x), by applying formulas in 
lH pp. 100-105]. 

Now we examine F(4 i) (x) . This function is easy to examine, because when g = exp(— Tr^ri), 
where n E Q+, the values of an are called singular moduli and can always be calculated in 
terms of algebraic numbers (HI p. 214]. Their values have been extensively tabulated. For 
example, many explicit evaluations of an can be found in [4, pp. 281-306]. Since a has an 
argument of = e^^'^l^ , we write a = ae4/x2. When x = 8, then a = ai = 1/2, and 
therefore 

- - ^ ^ / 



-OO m=— OO 




Similarly, when x = 4/^7, we have [i4i p. 284] a = ctzs = ( - 1)^ (2^2 - V7) 
m p. 284]. Thus we obtain 

V V (-1)"^" _ ^ 1 A + (v^-l)v/272W7 ^ 

(4m)2 + 7(4rz -fip 4v/T4 ^ 1^ 1 - (v/2 - 1)^2^2 - v^^ ' 



n=— OO r?i=— OO 



There are many similar identities that follow from results in [4], but they often tend to be 
very complicated. The majority of the identities contain algebraic numbers involving nested 
radicals. 

We conclude this section by proving a pair of formulas for -F(5,i) (x) . By (13.51) . this requires 
calculating the parameter fi. In principle, these calculations are straight-forward exercises. If 
the values of both R{—q) and R{q^) are known, then calculating n = R{—q)B?{q^) is trivial. 
If only one of the values is known, then fi can be calculated by solving (13.131) . This second 
type of calculation requires solving a cubic equation. If neither value is known, then we can 
use (|3.15l) to calculate ji. In practice, we have only been able to identify two instances where 
11 is reasonably simple. 

We begin by setting x = 1 in (1331) . By O pp. 57-58], 

^ = R{-e~^)R\e-^^) = ^(3 - V5)(7 -f 3^5) (^4 -f 2^5 - ^10(5+^5)^ . 
Substituting this last result into (13. 5|) and simplifying with Mathemat ica leads to 



Z E * ' 



ra=— OO ?Ti=— OO 

TT 



log (^-19 + 9^5 - 3^85 - 38^5^ + ^ log (yb - l) . 



(3.20) 



5^5 

Now we examine the more difficult case when x = 5. This choice of x leads to (11.11) 
quoted in our Introduction. We calculate /x using (|3.15l) and the values of class invariants 
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Gn, which are defined by 

where q = exp(— vryn). Hence, using the fact that Gn = Gi/n and the value of G25 
Bl p. 190], we find that 



1 - 4/i - ^5 26"^ 20^5 2 \^ 2 

Therefore jj, is given by 

/i = 1 (^3 + Ts) (^-4 + 2v^ - ^2(25-11^5)^ . 

Substituting this last result into (13.51) . and then simplifying nested radicals, we complete the 
proof of (fLTl) . 

4. Simplification for higher values 

In this section, we evaluate F(a,b){x) when a > 6. In order to simplify the calculations, 
we restrict our attention to cases where a E 2Z. In these cases we can apply elementary 
properties of Jacobian elliptic functions. Let us briefly recall that the elliptic functions sn(u), 
cn(M), and dn(M) are doubly -periodic, meromorphic functions, which depend implicitly on 
a parameter a = k"^, where k is called the elliptic modulus. Their periods are integral 
multiples of K and iK', where K and K' are complete elliptic integrals of the first kind 
associated with the moduli k and k' = a/1 — k"^, respectively. For us, the representations in 
terms of hypergeometric functions 2F1 E p. 102] 

A-:=^.F.(ii;l;a.). K':=^,F,(llia-a 

are employed in the sequel. There is a well-known inverse relation between the modulus, 
and the elliptic nome q given by jS] p. 102] 

tt = 1 jT^, g = exp -TT — ^ 

Theorem 4.1. Suppose that a G (0,1), and let 

(4.1) 

Assume that a G {2, 4, 6 . . . }, 6 G Z, and (a, h) = 1. Then 

a-l 

ax 



3=0 
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Proof. We have already established that (12.21 ) is true if a; G (0, oo) . The identity remains valid 
if we let 1 /a; — > 1/x + i. This substitution has the effect of sending q — )■ —q. Combining the 
two identities and performing a great deal of simplification leads to 

^-^ ^—^ [xmY + [an + 6)^ 



-12,+m, /^^ (l-^='+V"'+')(l-'^-''-'? 



2j-l 2m+l^ 



j=0 \m=0 ^ ^ ^ V ' 'i / 

Taking note of (14.11) . and then recalling the product representation for dn(M) [1 , p. 918], we 
find that the last expression transforms into 



oo 



f ^ < ^ ( T-m |2 _|_ ( r\nn -J- 



n=— oo m 



xmy + (an + 6)" 



j=0 

If we assume that a is even and 6 is odd, then the left-hand side of the identity equals 
2F(a,fe)(x). We can recover (14.21) by noting that J2j c^^^^-'^^^^ = whenever (a, 6) = 1. □ 

In order to provide an application of (14.21) . we evaluate F(io,i)(l) explicitly. While we re- 
strict our attention to this single example, the method we describe extends to many additional 
values of Fi^a,b)ix). Let us recall that dn(n) has real period 2K. If we use the symmetries 
[Hip. 500] 

dn{2K -u) = dn{u), dn{K - u) = , (4.3) 

dn(u) 

then (14.21) reduces to an expression involving [^] elliptic functions. When (a, b) = (10, 1), 
we have 




(4.4) 



-log(l - a). 



4 

By (|4.1I) . it is possible to calculate a whenever G and x > 0. It just remains to 
compute the values of the elliptic functions. 

Notice that dn(r_ft'/ s) is an algebraic function of a if (r, s) G Z^. This is a consequence of 
the fact that elliptic functions obey addition formulas [1 , p. 574]. Perhaps the easiest method 
for calculating values such as dn(f('/10) and dn(3i^'/10) is to generate polynomials (but not 
necessarily minimal ones) which they satisfy, by iterating the duplication formula for dn(z) 
m p. 574]. Let us recall that 

dn(2z) = fidn{z)), (4.5) 
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where 



X 



— [x^ — 1)(1 + -{x^ — 1)) 



For brevity, we use the shorthand notation 

di := dn ( ^ 



Using (14.31) and (14.51) . we can easily show that 



d2j = f{dj), 
dj = d2(}-j, 



_ VI - a 
"10- j - — -1 , 



d^ = \/\ — a. 

As a consequence of the elementary properties above, it is easy to deduce that 



- VT^ = 0, 



/(4)/(/(/(^^3)))- vr^ = o. 



(4.7) 



For instance, notice that f{di)f{f{f{di))) = d2f{f{d2)) = ^2/(^4) = rf2C?8 = Vl — It 
follows immediately that dn(i^'/10) and dn(3-ft'/10) are conjugate zeros of a function which 
is rational in Q(-\/l — a). It is easy to extract a polynomial which di and ^3 satisfy, by 
considering only the numerator of f(x)f{f{f{x))) — \/l — a. 

Now we can finish the computation of F(io,i)(l). Equation (14.11 ) shows that x = 1 when 
a = 1/2. As a result, (|4.4I) becomes 



n=— 00 m=— 00 



3j 

(4.8) 



5 -2v/5 , 
^ log 2, 



where 



rfi = dn(i\:/10) ^0.9915..., 



rfs = dn(3ir/10) ^ 0.9309 
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If we use Mathematica to expand (14.71) . then it is easy to see that di and are conjugate 
zeros of the irreducible polynomial 

=1 + 32X^ - 1152X^ + 14528X^ - 103328X^ + 445056X^° - 747008X^2 

- 5859584X1^ + 67132864X1^ - 404289024X1^ + 1770485760X^0 

- 6097568768X^2 ^ 17124502016X2^ - 40180561920X^6 + 80299532288X^8 

- 138787278848X^0 + 209592829440X^^2 _ 277574557696X3^ 
+ 321198129152X^6 - 321444495360X3® + 273992032256X^0 

- 195122200576X^2 ^ 113311088640X^^ - 51748995072X^6 
+ 17186013184X^8 - 3000107008X^0 _ 764936192X^2 

+ 911474688X^^ - 423231488X^6 ^ ii9013376X^8 _ 18874368X60 

+ 1048576X62 + 65536X6^ 

The calculation is essentially complete, but we provide a few additional comments. Despite 
the fact that Mathematica could not solve this equation directly, it is possible to express 
di and ds in terms of radicals, as we demonstrate below. It is unfortunate that the formulas 
are prohibitively complicated. 

We conclude by briefly describing how to recover explicit formulas for di and ^3. First 
notice that if f(x) = y, with f{x) defined in (14.61) , then we can express x in terms of y by 
solving quadratic equations. Since dn(2K/5) = /(/(rfi)) and dn(6-ft'/5) = /(/(ds)), it 
is sufficient to reduce dn(2_ft'/5) and dn(6i^/5) to radicals. There are several methods to 
accomplish this calculation. The simplest approach is to generate their minimal polynomials 
by repeated applications of the duplication formula for dn(z). It is then possible to verify 
the formulas 



In this example we have assumed that a = 1/2. It is still possible, albeit significantly 
more difficult, to evaluate these elliptic functions for certain other values of a. 



We have shown how to prove many explicit formulas for F(^a,b) (x) . The most obvious 
extension of this research is to examine cases where a > 5 is an odd integer. Notice that 
Theorem |4T| does not apply to those values. It should also be interesting to attempt to apply 
our techniques to the class of sums studied by Zucker and McPhedran in lfT6l . They gave 
closed form evaluations for many values of 




5. Conclusion 



00 



00 



,2s 




in terms of Dirichlet L-series. 
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